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The great east Japan earthquake and subsequent tsunamis caused Fukushima Dai-ichi Nuclear Power Plant 
(NPP) accident. National Institute of Radiological Sciences (NIRS) developed the external dose estimation 
system for Fukushima residents. The system is being used in the Fukushima health management survey. The 
doses can be obtained by superimposing the behavior data of the residents on the dose rate maps. For 
grasping the doses, 18 evacuation patterns of the residents were assumed by considering the actual 
evacuation information before using the survey data. The doses of the residents from the deliberate 
evacuation area were relatively higher than those from the area within 20 km radius. The estimated doses 
varied from around 1 to 6 mSv for the residents evacuated from the representative places in the deliberate 
evacuation area. The maximum dose in 18 evacuation patterns was estimated to be 19 mSv. 

The great east Japan earthquake occurred at 14:46, Japan time on 11 March 2011 in the Tohoku District, 
northeastern part of Japan. After the quake, big tsunamis rushed to the coast of the Tohoku district including 
Fukushima prefecture. The tsunamis caused electric failure of Fukushima Dai-ichi Nuclear Power Plant 
(NPP). The cooling systems of nuclear reactors stopped, and the hydrogen explosions occurred. Subsequently, 
radioactive fission products leaked out to the environment^'^. 

The interest in radiation doses due to the accident has been increasing among the Fukushima residents. The 
importance to estimate exposure doses of the residents has been recognized among the involved organizations. 
National Institute of Radiological Sciences (NIRS) started to develop the external dose estimation system for 
Fukushima residents at the end of March, 2011. On the other hand, Fukushima prefectural government and 
Fukushima Medical University decided to perform the health management survey for all Fukushima residents 
(about two million people) at the end of May, 201 P. External dose was considered to be one of the necessary data 
for the health management, and then the NIRS external dose estimation system was adopted in the survey. 

This paper describes the NIRS external dose estimation system and the calculation results for 18 evacuation 
patterns assumed by NIRS. 

Results 

We developed the NIRS external dose estimation system for the Fukushima residents to estimate the external 
effective doses for the first four months after the Fukushima Dai-ichi nuclear power plant accident. This system 
has been adopted in the Fukushima Health Management Survey, and the estimated result is sequentially provided 
to the individual Fukushima residents. 

The estimated external doses of the adult residents in 18 evacuation patterns were shown in Table 1. The results 
show that the doses of the residents evacuated from the representative places in the deliberate evacuation area 
(typically 1-6 mSv) were relatively higher than those from the area within 20 km radius (around 1 mSv). The 
maximum dose in 18 evacuation patterns was estimated to 19 mSv. 

Discussion 

The NIRS external dose estimation system for the Fukushima residents has been developed based on the 
simulation with the System for Prediction of Environmental Emergency Dose Information (SPEEDI)^ and 
monitoring data reported officially by the government. The uncertainties in this dose estimation depend mainly 
on not only those of the dose rate map but also the residents' memory on the behavior. 
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Table 1 Estimated external effective doses 


|mSv) for adult of 1 8 patterns 






Estimated external effective doses for adult (mSv) 




From the area within 20 km 


From the deliberate evacuation 


area 


Pattern No. 


Representative place 


Max. dose place Pattern No. Representative place 


k k 1 1 
Max. dose place 
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0.76 


0.76 


13 4.8 


13 
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0.94 


1.1 


14 0.84 


7.2 
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0.22 


0.64 


15 5.5 


17 
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0.23 


1.6 


16 6.2 


19 
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0.30 


0.30 


17 


1.8 
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0.25 


0.25 


18 


3.0 
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2.0 


2.3 
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0.77 
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1.5 






10 




0.18 






11 




0.74 






12 




0.71 







With regard to the dose rate map, since maps for the first 3 days are 
based on the SPEED I simulation instead of the monitoring data, the 
uncertainty is expected to be higher than those for other days. We 
have no alternatives for first 3 days at present because of the scarce 
monitoring data. The calculated results of SPEEDI greatly depend on 
the source term emitted from the nuclear reactors, which is expressed 
as a function of time. In the current situation, we have used the 
available and reliable data of the source term reported by the gov- 
ernment. However, when the accuracy improvement of the simu- 
lation is done with further analysis on the source term, recalculations 
with the dose rate maps base on the newly- simulated results would be 
necessary. On the other hand, the contributions of exposures in early 
stage for total estimated doses are not so much, because most resi- 
dents near the NPP have evacuated early after the accident and the 
increases of the dose by radionuclides in the plume were limited for a 
few hours as shown in Fig. 1. Actually the dose by radionuclides 
deposited on the ground etc., which is relatively plateau for the four 
months from the accident, is more important. For example, even at 



the Minami-soma city office where the increase of dose was observed 
earlier than other areas as shown in Fig. 1, the contribution of the first 
3 days to the total dose of four months was estimated to be about 5% 
by the NIRS system and about 10% when using the monitoring data. 
This fact indicated that the refinement of the simulation for the first 3 
days would not so affect to the total dose. 

In the NIRS system, the conversion coefficients from ambient dose 
equivalents to effective doses for isotropic irradiation were adopted. 
ICRP does not provide the conversion coefficients for the special 
condition of the irradiation. Saito et al.^ calculated the conversion 
coefficients of effective dose from environmental gamma ray sources: 
cloud and ground shines. The coefficients of isotropic irradiation 
(ISO), lateral irradiation from the right side to the left side of the 
body (RLAT) and lateral irradiation from the left side to the right side 
of the body (LLAT) are close to that assumed gamma ray source on 
the ground according to ICRP 74^ and Saito el al., and the coefficient 
of ISO is the most conservative. Table 2 shows the conversion coeffi- 
cients from ambient dose equivalent to effective dose for adult for 




□ Fukushima city (Momijiyama) 
Fukushima city (Oyama-chou) 

-•€>- litate village 

-■A- Tamura city 

-"©- Kawauchi village 
Iwaki city 

-•€>- Minamisouma city 
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Figure 1 | Monitoring data reported by MEXT and Fukushima prefectural government^'^'^. 
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Table 2 | Conversion coefficients from ambient dose equivalent to 
effective dose for adult for various radionuclides 

Nudide Xe-133 Te-129m Te-132 1-131 1-132 Cs-134 Cs-137 



Saitoetal.^ 
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0.61 


0.51 
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(Ground) 
















Saitoetal.^ 


0.34 


0.57 


0.46 
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0.57 


0.57 


(Air) 
















Present 


0.44 


0.58 


0.48 


0.53 


0.59 


0.58 


0.57 



various radionuclides considered in this study calculated from ICRP 
74 (adopted in this study) and Saito el al. This result indicated that 
conversion coefficient of 0.6 used in the NIRS system is reasonable. 

Methods 

Flow of the dose calculation. In the NIRS system, the external effective dose between 
March 12 and July 11, 2011 is estimated by superimposing the individual behavior 
data of each day on the daily dose rate map of such day. The data flow in the external 
dose estimation system is shown in Fig. 2. The flow is composed of three parts: 1) 
behavior data of residents, 2) construction of daily dose rate maps, and 3) calculation 
of effective dose. Each part is described below. This system utilizes two kinds of the 
fundamental data: daily dose rate map and behavior data of the resident. The 
background dose due to natural radiation sources was not included in the dose rate 
maps, because of the radiation protection purposes. The body size as a function of age 
is also considered because the external effective dose depends on body size even in the 
same radiation field. The dose rate maps have the data of effective dose for adults, so 
body size correction is considered for children less than or equal to 15 years old by 
using age coefficients. 

Data of behaviors of the residents. Fukushima prefectural government has been 
conducting the health management survey for the Fukushima residents to support 
their health conditions affected by the accidents in cooperation with Fukushima 
Medical University^ The survey consists of two kinds of investigations, which are 
called the basic survey and the detailed survey, respectively. In the basic survey. 



questionnaires regarding behavior of a resident for four months after the accident 
were sent by postal mail to the all residents in Fukushima prefecture, about 2 million 
in the population. Based on such data, the behavior data of the resident for each day 
can be obtained as follows; 

• Place to stay, /, 

• Time of stay in the building at /, ti^j, 

• Time of stay outside the building at /, toutj, 

• Type of the building at /, kj, 

• Time of the move from z to / + i, tmova-^i+h (if any) 

Construction of daily dose rate maps. The dose rate maps consisting of a lot of 
divisions were made based on the simulated and measured data in a time series 
between 12 March and 11 July in 2011. The size of each division in the map is 
corresponding to approximately 2 km X 2 km, which is based on the size divided the 
second mesh (latitude distance of 7.5 min. by longitude distance of 5 min) defined by 
Geospatial Information Authority of Japan. 

The daily dose rate maps were constructed based on two kinds of data. One is the 
hourly effective dose rate maps simulated by SPEEDI with the source term calculated 
by the MELCOR code^° by the Nuclear and Industrial Safety Agency (NISA), which 
was used between March 12 and 14, 201 r\ These data are an alternative to the 
monitoring data, because the number of measurement point was not sufficient to 
construct the dose rate maps in the period. The hourly effective dose rate maps 
outputted from the SPEEDI simulations, were averaged over a day, because the dose 
rate maps used in our system were daily maps. Also, its area was divided to those of 
1 km X 1 km. Therefore, the dose rate maps generated by SPEEDI were recon- 
structed by dividing in proportion to the area size of our system (approximately 2 km 
X 2 km) with commercially-available mapping software. As a result, the effective 
dose (for adult) rate maps with d{h,m), where h and m are the date and the division 
number, respectively, were obtained between 12 and 14 March, 201 1. Unfortunately, 
their area was limited in 98 km X 98 km, which is painted in green in Fig. 3, 
depending on the SPEEDI ability. At present, the doses are not added outside of the 
area. 

Another data is monitoring data released by the Ministry of Education, Culture, 
Sports, Science and Technology (MEXT), which was used between March 15 and July 
11, 2011. The monitoring data scattered in the map were converted to the 
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Figure 2 | Data flow in the NIRS external dose estimation system. 
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Figure 3 | Areas of dose rate maps used in the NIRS external dose estimation system. Basic maps were created by using ArcGIS of ESRI. 



spatially-continuous data by using Natural Neighbor method and then the daily dose 
rate in each division of approximately 2 km X 2 km was obtained by averaging the 
values in such division. Since the monitoring data maps by MEXT could not cover the 
area of a part of seaside and the boundary between Fukushima and Niigata prefectures 
(pink- colored areas in Fig. 3), the values of their neighbor on the right or left were 
used alternatively. 

Monitoring data on 15 March, 2011 were not sufficient to construct the dose rate 
map. On the other hand, this SPEED! results were significantly different from the 
monitoring data. This was expected to derive from the highest time variation of the 



dose rate in Fig. 1. Therefore, the dose rate map on 16 March was alternatively used 
for that on 15 March. Table 3 shows the comparison of effective doses on 15 and 16 
March, 2011 based on the monitoring data at the points shown in Fig. 1. The dif- 
ferences were between — 1 1 and + 182 |iSv (effective dose), so the alternative does not 
lead to a significant underestimation. In fact, it leads to an overestimation at most 
points. 

In the NIRS system, the background dose from natural radionuclides, etc. is not 
included in the estimated dose. Therefore, the background dose was estimated from 
the reported data before the accident. Fukushima prefectural government has 



Table 3 Comparison of effective dose 


(|iSv) on 1 5 and 1 6 March, 201 1 based on the monitoring data 




Effective dose (|iSv) 


Measurement point 


15 March, 2011 (A) 


16 March, 2011 (B) 


Difference (B-A) 


Fukushima city (Momijiyama) 


70 


230 


160 


Fukusima city (Oyama-chou) 


100 


265 


165 


litate village 


172 


354 


182 


Tamura city 


3 


5 


2 


Minamisouma city 


43 


58 


14 


Kawauchi village 


31 


23 


-7 


Iv/aki city 


65 


55 


-11 
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Figure 4 | Ratios of conversion coefficients for each age group to adult. 



reported the averaged air kerma rates per month at 23 points in the prefecture (1 in 
Kouno town, 4 in Naraha town, 5 in Tomioka town, 5 in Okuma town, 4 in Futaba 
town and 4 in Namie town) as 33-54 |j,Gy/h^^. The other measured data showed 
0.04 |iGy/h in Fukushima city, 0.04-0.06 |aGy/h in Koriyama city, 0.04-0.05 |aGy/h 
in Shirakawa city, 0.04-0.05 |aGy/h in Aizuwakamatsu city, 0.02-0.04 |aGy/h in 
Minami-aizu town, 0.05 |iGy/h in Minami-soma town and 0.05-0.06 |aGy/h in Iwaki 
city^. As the background dose rate, 0.04 |iGy/h as the median value among them was 
used for the background subtractions. Actually, the background effective dose rate 
was set to 0.03 |aSv/h by applying the conversion coefficient from air kerma to 
effective dose. 



The monitoring data are measured in terms of ambient dose equivalent H*(10), 
which is an operational guantity. On the other hand effective dose E is protection 
quantity and used for risk assessment of human health^^'^'*. Then the conversion 

- were used for 



factors from ambient dose equivalent to effective dose - , , 
^ H*(10) 

calculations in the system. According to the previous study^^, the main radionuclides 

emitted from the Fukushima Dai-ichi NPP were Xe- 133, Te- 132 and I- 132 at the early 

stage of the accident, and Cs-134 and Cs-137 afterwards. Adding 1-131 and Te-129m 

to these radionuclides, the conversion coefficients of effective dose per ambient dose 

equivalent for each radionuclide were calculated in the following steps; 

1) Calculate conversion factors of ( ^ j for isotropic irradiation (ISO) and 
H*(10)\ 

■ j for the energy of the i-th gamma ray emitted from a radionuclide, 

and {~~^^~^ ' according to the ICRP publication 74 with the logTog inter- 
polation method. Here, (]) is the fluence of gamma ray. 



2) Calculate conversion factors of 



H*(10) 



for each nuclide. 



and 



for isotropic irradiation (ISO) and 
' , by the tollowmg 



equations: 



E.((f). 



xR, 



T..R. 

^ I f H'jlO) 



xRi 



(1) 



(2) 



where Ri is the emission ratio of the i-th gamma ray emitted from a radionuclide. 



3) Divide ( - J by 

nuclide 



H*(10) 




Figure 5 | Assumed 18 evacuation patterns of the residents. Basic maps were created by using ArcGIS of ESRI. 
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The calculated — — r for each nuclide is shown in Table 2. In our estimation, the 
H*(10) £ 

value of 0.6 was adopted as — - — - for adult, rounding the value of 0.59 which is the 
H*(10) 

maximum among the radionuclides. Actually, since there are differences of tissue 
weighting factors of effective dose between ICRP Publication 60 and 103 esti- 
mated effective doses is not the same exactly depending on the publication. However, 
according to the reference^*^, these values are similar for isotropic irradiation. Hence 
the correction between them was not considered in the system. 

Finally, the effective dose (for adult) rate maps with d {h,m), can be obtained in a 
time series from 15 March to 1 1 July, 201 1 in all parts of Fukushima prefecture and a 
part of neighboring four prefectures (Miyagi, Yamagata, Tochigi, and Ibaraki) 
painted in pink and blue in Fig. 3. 

Calculation of effective dose. The external effective doses were calculated for three 
different situations: staying indoors/outdoors and moving from one place to another. 

When one stays indoors, the dose reduction should be considered because build- 
ings have shielding effects against radiation exposures depending on their materials 
and thickness of the walls. In the TECDOC 225^^ by the International Atomic Energy 
Agency (IAEA), representative reduction factors for cloud, r^-, and ground sources, 
are shown. From 12 to 14 March, the radio nuclides in the plume have contributed to 
the dose rates in the environment. On the contrary, the radionuclides on the ground 
were major sources of the exposures because of the rain or snow falls on 15 March in 
Fukushima. In the external dose calculation system, the reduction factors for cloud 
source in TECDOC 225 were adopted for the calculation from 12 to 14 March. The 
values used in this system were 0.9, 0.6, and 0.2 for one or two story wood frame 
house, one or two story concrete house or building, and three or more story concrete 
building, respectively. For the period of 15 March and after, the factors for deposited 
radioactivity were applied in the calculations. The values were 0.4, 0.2 and 0.1 for one 
or two story wood frame house, one or two story concrete house or building, and three 
or more story concrete building, respectively. 

Dose rates during a move may change depending on the location. The dose exposed 
during the move is simply calculated as the product of mean effective dose rate at the 
start and end points multiplying the hours of the move. 

Finally, when the mesh number including the place, j is m(j), and dose reduction 
factor correspond to the type of the building, kj is rc/g(kj), dose for stay on the date, h, 
Estay,h and dose for the move on the date, h, E^ovch express as the following equations, 
respectively. 

EstyM = H^Mj)) X [tinj X r,/g{k,) + tout,j] } (3) 
Emove,h = 2^. < X tmove,M+l > (4j 

The effective dose from March 12 to July 11, 2011, £, can be obtained with the 
following equation. 

E= ^ ^ {Estay,h~\' Efnove,}i) (5) 
h 

The conversion coefficients from ambient dose equivalent to effective dose El 
H*(10) are higher as the age is lower, because the shielding effect increases with the 
volumes of body stem^*^. When the maximum conversion coefficient was used for a 
conservative estimation, conversion coefficient. Cage (age coefficient) was expressed 
as the following linear function of age, y, between 1 and 15 years old (Fig. 4). 

Cage(y)= -0.0144x7+1.27 (6) 

The effective dose for children can be calculated by multiplying E by Cage- 

Eighteen patterns of the evacuations. For roughly grasping and estimating the doses, 
18 evacuation patterns of the residents from the area within 20 km radius of the NPP, 
and the deliberate evacuation area^ were assumed by considering the actual 
evacuation information before using the survey data (Fig. 5). 

Since the traffic conditions were worse compared to those in ordinary situation, the 
moving speeds of evacuations were assumed to be 10 km/h for the moving inside the 
prefecture, and 30 km/h for those from Fukushima to other prefectures. Until the end 
of March, the ordinary traffic conditions were considered. For the time of activities 
outside the houses or buildings, it is assumed that the residents stayed for 4 hours a 
day before the evacuations, and 2 hours a day after the evacuations. 

The two places where the residents started to evacuate were assumed; the repres- 
entative place such as town office, and/or the places where the doses were highest in 
the town until when they evacuated. 
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